%dXt = sin(5t)dW 
clear all
% index J for X
alphaK=15; alphas=[zeros(1,alphaK);eye(alphaK)];
% time T
T = 1/2;
% time mesh 
num=50; tspan = linspace(0,T,num+1);
%sigma
s=mihandle;
s.sigma=@(s) sin(5*s);
%options = odeset('AbsTol',1e-12,'RelTol',1e-6);
options = [];
row=size(alphas,1); X0=zeros(row,1);
% solver of X
disp('solving X')
tic
[time,X] = ode45(@randomsg,tspan,X0,options,T,alphaK,s);
toc
varitest(time,X.*X,T,'single','sin5t',num)